NASA TN D-6930 


lOJSi- 30X&0 


NASA TECHNICAL NOTE 



NASA TN D-6930 


CASE F* LE 

COP-Y 


OPTIMUM RUNWAY ORIENTATION 
RELATIVE TO CROSSWINDS 


by L. W. Falls and S. C. Brown 
George C. Marshall Space Flight Center 


Marshall Space Flight F enter, Ala. 35812 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • SEPTEMBER 1972 



TECHNICAL REPORT STANDARD TITLE PAGE 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

MSA TN D-6930 

4. TITLE AND SUBTITLE 

Optimum Runway Orientation Relative to Crosswinds 


7. AUTHOR (S) 

L.W. Falls and S. C. Brown 

9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


12. SPONSORING AGENCY NAME AND ADDRESS 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


3. RECIPIENT'S CATALOG NO. 


5. REPORT DATE 

September 1972 


6. PERFORMING ORGANIZATION CODE 


8. PERFORMING ORGANIZATION REPORT # 
10. WORK UNIT NO. 

1 1 CONTRACT OR GRANT NO. 

13. TYPE OF REPORT & PERIOD COVERED 


[ 14. SPONSORING AGENCY CODE 


15. SUPPLEMENTARY NOTES 

Prepared by Aero-Astrodynamics Laboratory, 

Science and Engineering 

16. ABSTRACT 

Specific magnitudes of crosswinds may exist that could be constraints to the success of an 
aircraft mission such as the landing of the proposed space shuttle. A method is required to 
determine the orientation or azimuth of the proposed runway which will minimize the proba- 
bility of certain critical crosswinds. Two procedures for obtaining the optimum runway 
orientation relative to minimizing a specified crosswind speed are described and illustrated 
with examples. The empirical procedure requires only hand calculations on an ordinary 
wind rose. The theoretical method utilizes wind statistics computed after the bivariate normal 
elliptical distribution is applied to a data sample of component winds. This method requires 
only the assumption that the wind components are bivariate normally distributed. This 
assumption seems to be reasonable. Studies are currently in progress for testing wind 
components for bivariate normality for various stations. The close agreement between 
the theoretical and empirical results for the example chosen substantiates the bivariate 
normal assumption. 



19. SECURITY CLASSIF. (of thU report! 

Unclassified 


20. SECURITY CLASS1F. (of thU page) 

Unclassified 


21. NO. OF PAGES I 22. PRICE 


MSFC - Form 3292 (May 19 


* For sale by the National Technical Information Service. Springfield, Virginia 22151 











TABLE OF CONTENTS 


Page 

INTRODUCTION 1 

EMPIRICAL METHOD 1 

THEORETICAL METHOD 6 

CONCLUSIONS 9 

REFERENCES 11 

BIBLIOGRAPHY 11 

APPENDIX. COMPUTER PROGRAM FOR OPTIMUM RUNWAY 

ORIENTATION RELATIVE TO CROSSWINDS 12 

ACKNOWLEDGMENT 

The authors wish to acknowledge the assistance of Mr. Raymond Smith 
of the Computation Laboratory, MSFC, for his work on the Univac 1108 
Computer Program. 


iii 


TECHNICAL NOTE D- 


OPTIMUM RUNWAY ORIENTATION 
RELATIVE TO CROSSWINDS 

INTRODUCTION 


Runway orientation is undoubtedly influenced by a number of factors — 
perhaps winds, terrain features, population interference, etc. In some 
cases, the frequency of occurrence of crosswind components of some signi- 
ficant speed may have received insufficient consideration. If, for example, 
a runway for the space shuttle vehicles is being planned, it may be prudent 
to consider the optimum runway orientation to minimize crosswind com- 
ponents of, for example, 20 knots. Aligning the runway with the prevailing 
wind will not ensure that crosswinds of this magnitude will be minimized. 

In fact, two common synoptic situations (one producing light easterly winds, 
and the other causing strong northerly winds) might exist in such a relation- 
ship that a runway oriented with the prevailing wind might be the least use - 
ful to an aircraft constrained by crosswind components = 20 knots . Two 
methods (one empirical, the other theoretical) of determining the optimum 
runway orientation to minimize critical crosswind component speeds are 
described below. Both methods gave identical results for the Cape Kennedy, 
Florida, winds shown in Table 1. 

EMPIRICAL METHOD 

The following paragraphs outline a short procedure (one requiring 
only a desk calculator and a wind rose) for determining the best runway 
orientation relative to some specified crosswind component. 

From the ordinary wind rose (Table 1) , the percentage frequency 
of a number of wind speeds and directions can be obtained. For this pro- 
cedure, the percentage frequency or number of cases in each class interval 
or "box" is assumed to be located at the class mark. 

It is apparent that, if two of the reported wind directions (opposites) 
are chosen as a runway orientation, the crosswind component contributed by 
each box can be obtained from the product of the wind speed and the sine 


TABLE 1. ANNUAL WIND ROSE AND STATISTICAL SUMMARY FOR CAPE KENNEDY 
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TABLE 1. (Concluded) 
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the major axis of the distribution. 




of the angle between the wind direction and the runway orientation. Since 
wind directions are usually given to 16 points, the angles between the wind 
direction and the runway will be 22. 5, 45, 67. 5, and 90 deg. For example, 
if the runway orientation is chosen as EW, then all ENE, ESE, WNW, and 
WSW winds will be 22. 5 deg off the runway. Likewise, NE, SE, NW, and 
SW will be 45 deg off, NNE, SSE, NNW, and SSW will be 67. 5 deg off, and N 
and S will be 90 deg off. From Figure 1 it is apparent that each wind direction 
and wind speed will make a contribution to the total percent frequency of cross- 
wind components. Since several wind directions in the same wind speed category 



Figure 1. Runway crosswind component calculation. 
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will produce the same crosswind component, these boxes should be summed. 
For example, under the 10-19 category, sum NE, SE, NW, and SW (2180 + 
3841 + 2866 + 1548 = 10 4 35) to obtain the frequency of crosswind components 
of 10. 3 knots ( 14. 5 sin 45 deg = 10. 3) . This process should be continued until 
all boxes except E-W have been used to compute a crosswind component. 

In summary, the crosswind component computation procedure consists 
of the following steps: 

1. Compute all possible crosswind components — the product of each 
wind speed and the sine of 22. 5, 45, 67. 5, and 90 deg. 

2. Sum all boxes that contribute the same crosswind component for a 
specified runway orientation. Compute the frequency and percent frequency 
for each crosswind speed. 

3. Order the crosswind component speeds from the largest to the 
smallest and tabulate the percent frequency of occurrence opposite each cross- 
wind component. 

4. Form the cumulative percentage frequency (CPF) from the values 
tabulated in step 3. , starting with the highest wind speed. This CPF gives a 
description of the crosswind components for a single runway. 

5. Interpolate the CPF for the desired wind speed. This interpolated 
value gives the probability in percent of equaling or exceeding the specified 
crosswind component for that runway orientation. From the few cases 
examined, it appears that the interpolation should be made assuming a normal 
distribution. 

Of course, the procedure must be repeated for each pair of opposite 
wind directions (runway orientations) to determine the optimum runway 
orientation relative to a critical crosswind component. 

This procedure was applied to the Cape Kennedy annual wind rose 
shown in Table 1, except that hurricane associated winds > 50 knots were 
removed on the premise that landing operations would not be conducted during 
such periods. (Table 1 was prepared from hourly peak wind measurements 
made from September 1958 through June 1969. ) Results of the analysis 
( Fig. 2) , for which 20 knots was selected as the critical crosswind speed, 
indicate that the best runway orientation relative to speeds of this magnitude 
is about 150 to 330 deg true. 
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180 200 220 240 260 280 300 320 340 360 

Figure 2. Cape Kennedy percentage frequency of exceedance of peak 
crosswind components > 20 knots versus runway orientation. 

THEORETICAL METHOD 

For Cape Kennedy, the assumption has been made that wind components 
are bivariate normally distributed; i.e. , a vector wind data sample is resolved 
into wind components in a rectangular coordinate system and the bivariate nor- 
mal elliptical distribution is applied to the data sample of component winds. 

For example, let x t and x 2 be normally distributed variables (the wind compo- 
nents at Cape Kennedy) with parameters (£ t , c^) and (| 2 , o- 2 ) . i \ and £ 2 are the 
respective means, while cti and cr 2 are the respective standard deviations. Let 
p be the correlation coefficient, which is a measure of the dependence between 
Xj and x 2 . Now, the equation of the bivariate normal density function is 



Let a be any arbitrary angle in the rectangular coordinate system and 
let the sample estimates of £ lt £ 2 , cj, a 2 , and p be denoted by Xj, X 2 , S , 


S , and r , respectively. 
X 2 XjX 2 


From the statistics in the (x 1# x 2 ) space, the 


statistics for any rotation of the axes of the bivariate normal distribution 


through any arbitrary angle a may be computed as given below [1] . 
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Let (yi, y 2 )_ be the new space after rotation with orthogonal axes y t 

and Vo. LetYi, Y„ S , S , and r be the means, standard deviations, 
J 1 1 * y y y y 

12 12 

and correlation coefficient in the (y , y ) space. From Reference 1, 

Yj = X t cos a + X 2 sin a 

Y 2 = X 2 cos a - Xj sin a 

S = IS 2 cos 2 q/+S 2 sin 2 a; 

Yi \ x i x 2 


S = ( S 2 cos 2 q; + S 2 sin 2 o! 
y 2 \ x 2 x l 

S 

r = yiy?- 

yiy 2 S S ’ 

2 Yi y 2 

where S 2 , S 2 , and S are the variances and covariance in the (x,, x.) 

v V v V ' A ^ 

X 1 x 2 x l x 2 

space and S is the covariance in the (yi, y 2 ) space. 
yiy2 

In equations ( 2) , all statistics on the right side are available for Cape 
Kennedy or may be computed from existing data ( Table 1) . 

Equations (2) give the statistics in the (y u y 2 ) space (this is the space 
after rotation through any angle a) that defines the bivariate probability den- 
sity function (1). 

The existing wind statistics for Cape Kennedy ( Table 1) are for the 
reference angle a = 90 deg; i.e. , the (x lf x 2 ) space is for a = 90 deg. 

Let Ao; denote the desired increments for which runway orientation 
accuracy is required; e.g. , one may wish to minimize the probability of cross- 
winds with a runway orientation accuracy down to Aa = 10 deg. This means 
we must rotate the bivariate normal axes through every 10 deg beginning at 
a = 30 deg [the (x lt x 2 ) space]. It is only necessary to rotate the bivariate 
normal surface through 180 deg, since the distribution is symmetric in the 


l 2 


+ 2S cos a sin 
x l x 2 


\l/ 2 
in a J 


\ l/ 2 

— 2S cos a sin a \ 
x l x 2 / 


( 2 ) 
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other two quadrants. This process will result in 18 sets of statistics in the 
(Yi» Y2) space obtained from equations (2) . y t is the head wind component 
while y 2 is the crosswind component. Since we are concerned with minimizing 
the probability of crosswinds (y 2 ) only, we now examine the marginal distri- 
butions p(y 2 ) for the 18 orientations ( oi ). Since p(y lt y 2 ) is bivariate normal, 
18 marginal distributions p(y 2 ) must be univariate normal: 

P( Y 2 ) = [cr 2 ( 2tt) */*] _1 exp | [(y 2 - £ 2 ) /a 2 ] 2 | . (3) 

£ 2 and cr 2 are replaced by their sample estimates Y 2 and S . Now, let 

y 2 



(4) 


where y 2 is the critical crosswind of interest, z is a standard normal variable 
and the probability of its exceedance is easily calculated from tables of the 
standard normal integral. Since a right or left crosswind (y 2 ) is a constraint 
to an aircraft, the critical region ( exceedance region) for the normal distribu- 
tion is two-tailed; i. e. , we are interested in twice the probability of exceeding 
| y 2 1 . Let this probablity of exceedance or risk = R. 

A computer program is available for computing the statistics in the 

(y lt Y 2 ) space defined by equations (2), the integral of equation (3), and the 

risk R (see Appendix A) . The Cape Kennedy data given in Table 1 are used as 

an example. The input to the program is Xj = -1.775, X 2 = -0.210, 

S = 9. 230, S = 9. 658, r = -0. 145, and Ace = 10 deg. Let the crit- 
x t X 2 XjX 2 

ical crosswind component x 2 = 20 knots as in the empirical example in the 
preceding section. To obtain a valid comparison to the empirical probability, 
we will compute R = the probability of exceeding or equalling x 2 = 20 knots. 
Table 2 summarizes the runway orientation angles a, the standard normal 
variate z, and the risk R for Cape Kennedy. 

For this example, the risk R = 5. 63 percent of exceeding or equalling 
a 20-knot crosswind is a minimum occurring at the runway orientation 
a = 150 - 330 deg. This theoretical result verifies the empirical conclusion 
arrived at in the previous section. 
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TABLE 2. SUMMARY OF RUNWAY ORIENTATION (a ) VERSUS 
RISK (R) , CAPE KENNEDY. 


a 

(deg) 

z 

R 

(%) 

a 

(deg) 

z 

R 

(%) 

90 

1.57 

11.53 

180 

1. 82 

6.91 

100 

1.65 

9.93 

190 

1. 76 

7.79 

110 

1.72 

8.47 

200 

1.71 

8.77 

120 

1.80 

7.25- 

210 

1. 66 

9.76 

130 

1.86 

6.36 

220 

1.61 

10.69 

140 

1. 89 

5.82 

230 

1.58 

11.48 

150 

1 

1.91 

5.63 

240 

1. 55 

12.07 

160 

1.90 

5.78 

250 

1. 54 

12.45 

170 

h- 1 

• 

00 

Oi 

6.22 

260 

1.53 

12. 61 


The method described to determine the orientation of a runway which 
will minimize the probability of critical crosswinds is accurate and expedient. 
The procedure described may be used for any station. Only parameters esti- 
mated from the data are required as input. Consequently, many runways and 
locations may be examined rapidly. 

CONCLUSIONS 


Either the empirical or theoretical method described in this report 
may be used to determine an aircraft runway orientation that minimizes the 
probability of critical crosswinds. Again, it is emphasized that the wind 
components must be bivariate normally distributed to use the theoretical 
method. The Cape Kennedy wind component raw data were not available for 
the example used in this report. The bivariate normal assumption may or 
may not be a good assumption. However, the Cape Kennedy sample illustrates 
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the method. The agreement between the theoretical and empirical methods 
is very good. As shown by Figure 2 and Table 2, the minimum and maximum 
probabilities occur at the same azimuths for both methods. The minimum 
risk of 5. 63 percent for the theoretical method occurs at 150 deg. The mini- 
mum risk for the empirical approach also occurs at 150 deg as shown by 
Figure 2. The empirical risk is approximately 4. 7 percent. The maximum 
risks for the two methods also occur at approximately the same orientation. 

The maximum theoretical risk of 12. 61 percent occurs at 260 deg, while the 
maximum empirical risk of about 7. 6 percent occurs at 250 deg. One may 
view the differences in the theoretical and empirical probabilities as a measure 
of the departure of the data from normality. 

In practical applications, the following steps are suggested: 

1. Test the component wind samples for bivariate normality if these 
samples are available. See Reference 2 for bivariate normal goodness-of-fit 
tests. 1 


2. If the component winds are available and cannot be rejected as 
bivariate normal using the bivariate normal goodness-of-fit test, use the 
theoretical method since it is more expedient and easily programmed. 

3. If the component wind data samples are not available and there is 
doubt concerning the assumption of bivariate normality of the wind components, 
use the empirical method. 


1. H. L. Crutcher and L. W. Falls: 
published) . 


Multivariate Normality. 1971 (to be 
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APPENDIX. COMPUTER PROGRAM FOR OPTIMUM RUNWAY 
ORIENTATION RELATIVE TO CROSSWINDS 

Legend 


oi Runway orientation 


a (R) Orientation or angle at which the wind component statistics are 

computed from the data sample 


Xj Mean of the head-tail winds at orientation a (R) 


X 2 Mean of the crosswinds at orientation a (R) 


S Standard deviation of head-tail winds at ck (R) 

x i 


S Standard deviation of crosswinds at oi (R) 

x 2 


S Covariance at q;(R) 

x l x 2 

R Correlation coefficient at a (R) 

x l x 2 


Yj Mean of the head-tail winds at orientation a 


Y 2 Mean of the crosswinds at orientation a 


S 2 Variance of the head-tail winds at q; 

Yi 


S Standard deviation of the head-tail winds at a 

yi 


S 2 Variance of the crosswinds at a 

y 2 


S Standard deviation of the crosswinds at a 

Y2 


S Covariance at orientation a 

yiY2 
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R 

y iY2 


Correlation coefficient at a 


y 2 Critical crosswind 

z Standard normal variable 

F (y 2 ) Normal distribution function 

R Risk of exceeding | y 2 | 

INPUT 

SALP Starting a 

ALPR aR 

XB1 Xj 

XB2 X 2 

SX1 S 

X 1 

SX2 S 

~ x 2 

RX1X2 R 

x l x 2 

DALP Aa 

EALP End or Maximum a 

NV No. of Y2 

Y2 Table of Y2 

The first Y2 is taken and the program is utilized with the starting a. 
a is then increased by DALP and the program is utilized again. This is 
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repeated until EALP (maximum a) is reached. The next Y2 is taken and the 
process is repeated. This is continued until all of the Y2's have been used. 

A transfer is then made to STATEMENT to see if more data are to be read 
in. If not, the program terminates. Usually Y2 is the only data that vary. 

Note: The STATEMENT following STATEMENT 10 is an end-of-file 
check. If there are no more data the program terminates. This end-of-file 
test may not be compatible with other computers. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, June 2a, 1972 
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